The Millennium Galaxy Catalogue: The nearby supermassive black hole 
mass function 
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Abstract. We highlight the correlation between a galaxy's supermassive black hole mass and the 
Sersic-index of the host spheroid or bulge component. From our bulge-disk decompositions of 10 095 
galaxies, drawn from the Millennium Galaxy Catalogue, we construct the local {z < 0.18) mass function 
of supermassive black holes. We compare our results to those of McLure & Dunlop (2004) and conclude 
that the mass density of supermassive black holes may be marginally higher than previously supposed. 
Oh This increase is predominantly due to the inclusion of low mass and later-type bulges. More details will 

' be presented in a forthcoming paper. 

CO ' 1 Introduction 

<n : 

Supermassive black holes (10 9 — 10 6 M Q ) have now been identified through core velocity measurements for ~ 40 
nearby systems, see for example the compendium of Ferrarese & Ford (Hj- The formation mechanism of these 
^-j- ■ super massive black holes remains uncertain. Silk & Rees |2.'ij advocate formation directly through initial 

monolithic collapse whereas Croton et al. [2] advocate a gradual build-up through successive mergers. It is 
£ — , also unclear whether central black holes occur within all spheroidal, bulge and possibly nucleated systems. In 

Q\ ' particular, their existence (or the existence of intermediate- massive black holes), in dwarf systems is currently 

being pursued with no definitive result as yet, although see |2Z]. As well as their recent detection we are also 
now aware of the tight correlation between the black hole mass and the global velocity dispersion (a) of the 
bulge or spheroid which they inhabit |7|,|IH!- From the L — a relation [H] it also follows that there is a correlation 
"*^h \ between the black hole mass and the bulge- luminosity, albeit possibly non- linear |17j . 

Here we highlight a third correlation between the concentration of a galaxy bulge or spheroid and the 
supermassive black hole mass ^2]- This correlation is to be expected given the known relation between the 
luminosity, velocity dispersion and concentrations of nearby spheroid and bulge systems |26| . |13| . Which of 
4-> ■ the three relations (Mbh — o~, Mbh — L or Mbh — n) is the more fundamental is not yet clear. In practical 

terms the Sersic index is arguably the easier to measure as it a purely photometric quantity requiring neither 
. . ' kinematic observations (a) nor absolute photometric calibration (L). 

> ■ 

5h ' 2 Concentration and the Sersic index 

a ; 

The concentration of a galaxy's light is typically given as the ratio of two radii. For example the SDSS adopt 
the radius which encompasses 90% of the total flux divided by the radius which encompasses 50% of the total 
flux, i.e., 

(--(90/50) - r^- 

Alternate definitions are in common practice which use different ratios (e.g., £Q who use C 80 /2o f° r example). 
The concentration index can be related analytically to the Sersic index (see recent review by |15p of the projected 
light-profile ^3], where the Sersic intensity model ([21], |22])> is given by: 

J(r)=7 exp[-(p( 1 /")] ) 

where I(r) is the intensity at radius, r, I the central intensity, h the scale-length, and n the Sersic index or profile 
shape parameter. Note that n — 1 refers to the traditional exponential profile |§| and n = 4 the traditional de 
Vaucouleurs profile A value of n = 0.5 represents a Gaussian-shaped profile (see Fig. QJ. 



3 The SMBH-n correlation 



Fig.^left) shows the MsH — n relation while ¥\g\^right) shows the Mbh — o~ relation for the same systems. The 
linear fits are by ordinary least squares regression assuming a 20% uncertainty on n. Both distributions show a 
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Figure 1: Illustration of the projected Sersic profile with index varying from 0.5 to 10 as indicated. All profiles 
have the same half-light radius and surface brightness at the half-light radius but differing magnitudes. 



strong correlation with comparable (linear) Pearson, r, and (non- linear) Spearman, r s , rank-order coefficients. 
Work is underway to add the remaining dozen or so systems to this plot for which reliable black hole masses and 
distances are accurately known. Two further issues are worth mentioning. (1) The Sersic-index measurements 
are predominantly made in the R filter. While this is less susceptible to dust attenuation/distortion (and 
particularly for the spiral galaxy bulges), near-IR profiles would be more optimal as the Sersic index can be 
wavelength dependent (because of dust which is more centrally concentrated than the stars). (2) The sample 
is not entirely random but reflects those systems sufficiently nearby such that black hole mass constraints can 
be obtained. As we live in a specific region of the Universe (a loose group) this may bias our sample in some 
unforeseeable way. However, as the loose group environment appears to be the most common jS] any bias is 
likely to be small and subtle. The best linear fit to our data is: 

log M BH = (3.03 ± 0.48) log(n/3) + 7.82(±0.07). (1) 

This relation can be used to construct the SMBH mass function for any sample of galaxies for which the Sersic- 
index has been measured. As measuring the Sersic-index is much easier than the (bulge) velocity dispersion, 
which can be particularly problematic to measure for late- type systems due to disk/bulge contributions, one 
might expect the mass function derived via this method to be more robust. 




Figure 2: (left) supermassive black hole mass versus Sersic index of the bulge component or spheroid, and (right) 
supermassive black hole mass versus bulge or spheroid velocity dispersion (see section 3). Nonlinear fits are 
explored in Graham, et al (2005, in prep.) 

4 The MGC and bulge-disk decomposition 

The Millennium Galaxy Catalogue (MGC, |16| . 0]) has been discussed in two earlier talks at this meting (see 
Allen and Liske these proceedings). It constitutes a 37.5 sq degree imaging survey along the Northern Spring 
equatorial strip with redshift information for 96% of the galaxies for B < 20 mag. The sample has a median 
redshift of 0.12 and bulge-disk decomposition using GIM2D [21] has now been completed for all 10 095 systems 
(Allen et al. these proceedings). The survey has been used to constrain estimates of the nearby luminosity 
function ^1 as weu as measure the joint luminosity-surface brightness distribution Plans are afoot to expand 
the MGC program with additional ugriz data from VST and YJHK S data from UKIRT/ VISTA. The redshift 
survey will be extended to B = 22 mag with AAT/AAf2 and eventually to B = 24 mag with Gemini/WFMOS. 



5 The MGC-SMBH mass function 

To construct the SMBH mass function we need to combine the relation shown in Eqn. 1 with the luminosity 
function of galaxy bulges (see Liske et al. , these proceedings) . Our approach is to assign a black hole mass to 
each galaxy (based on the Sersic-index) and then adopt a weight for each galaxy dependent on the space-density 
of galaxies of that luminosity divided by the number of galaxies contributing to that calculation, i.e., 

W t = J>(M, »)/N(M, fj) = 1/Vmax- (2) 

Having established a list of weights, Wi, and masses, one can then simply sum the product of these distributions, 
thus: 

N(M BH )5\ogM BH = W(M BHi )i8\ogM B H. (3) 

i 

The result, shown in Fig. |3J includes genuine bulges but excludes pseudo-bulges (distinguished by their blue 
colours, see Liske et al. these proceedings). Also shown are recent estimates by based on M BB — o (dashed 
line) and M B h—L (solid line) estimates for a sample of ~ 9000 galaxies drawn from the Sloan Digital Sky Survey. 
The integral of this value leads to the total density of SMBHs which we tentatively find to be p = 3.4 x 10 5 
M q /iq 7 Mpc -3 . This is somewhat higher than the McLure & Dunlop values most likely due to our inclusion of 
later-type bulges. 



6 Conclusions 

We have highlighted the correlation between the mass of the super massive black hole and the Sersic index 
(i.e., the M BB — n relation). This is first presented in ^2] an d followed up in [T3J. The Sersic index of 
spheroids and bulges is straightforward to measure and we have now done so for 10 095 galaxies drawn from 
the Millennium Galaxy Catalogue. Using the M BB — n relation combined with our recent luminosity function 




Figure 3: The two lines show the SMBH mass functions derived by McLure & Dunlop (2004) using either the 
Mbh — & {dashed line) or the Mbh — L relation (solid line). Our results are show as data points with errorbars 
and appear to follow the Mbh — o results of McLure & Dunlop reasonably closely. 



estimate 0] we have now constructed the nearby mass function of super massive black holes. The shape of the 
distribution closely follows that found by JHj albeit with a slightly higher normalisation. We ascribe this to our 
method including all bulges whereas |19| exclusively focused on spheroids and early- type spiral systems. The 
result presented here will be presented in a more robust fashion in the forthcoming paper by Graham et al., in 
preparation. 
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